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ABSTRACT 

The incidence and properties of Active Galactic Nuclei (AGN) in the field, groups, and clusters 
can provide new information about how these objects are triggered and fueled, similar to how these 
environments have been employed to study galaxy evolution. We have obtained new XMM-Newton 
observations of seven X-ray selected groups and poor clusters with 0.02 < z < 0.06 for comparison with 
previous samples that mostly included rich clusters and optically-selected groups. Our final sample 
has ten groups and six clusters in this low-redshift range (split at a velocity dispersion of a = 500 
km s _1 ). We find that the X-ray selected AGN fraction increases from /_4 (Ljf > 1Q 41 ;Mr < Mj£+1) = 



0.047 



+0.023 
-0.016 



in clusters to 0.091 



+0.049 
-0.034 



for the groups (85% significance), or a factor of two, for AGN 
above an 0.3-8keV X-ray luminosity of 10 41 erg s _1 hosted by galaxies more luminous than The 
trend is similar, although less significant, for a lower-luminosity host threshold of Mr = —20 mag. For 
many of the groups in the sample we have also identified AGN via standard emission-line diagnostics 
and find that these AGN are nearly disjoint from the X-ray selected AGN. Because there are substantial 
differences in the morphological mix of galaxies between groups and clusters, we have also measured 
the AGN fraction for early-type galaxies alone to determine if the differences are directly due to 
environment, or indirectly due to the change in the morphological mix. We find that the AGN fraction 
in early-type galaxies is also lower in clusters fA,n>%s(Lx > 10 41 ;Mr < M£ + 1) = 0. 048^^019 
compared to 0.119±g;g^ for the groups (92% significance), a result consistent with the hypothesis 
that the change in AGN fraction is directly connected to environment. 

Subject headings: galaxies: active - galaxies: clusters: general - galaxies: evolution - X-rays: galaxies 
- X-rays: galaxies: clusters - X-rays: general 



1. INTRODUCTION 

While it is clear that the Active Galactic Nuclei are 
powered by accretion onto supermassive black holes 
(SMBHs), and that this accretion requires both a source 
of fuel and a trigger to remove angular momentum 
from the fuel, the origin of the fuel and the trigger 
mechanism(s) remain poorly understood. For lumi- 
nous QSOs, major mergers between gas-rich galaxies 
are considered the dominant mechanism and are the 
best, and perhaps only, candidate. Many studies 
have shown that a large fraction of these galaxies are 
morphologically disturbed, with close neighbors, tidal 
tails, multiple nu clei, or linked by luminous matter to 
other gal axies dGehrens. Fried. WetTinger, fc Wyckofil 
1984 iHutchings, Crampton . Campbell 11984 



Malkan. Marg on. fc Chananl 11984 iSmith et al.l 11986 
These results have motivated the hypothesis that 
lower-luminosity AGN in the local Universe are, like 

1 Current Address: Steward Observatory, University of 
Arizona, 933 N. Cherry Ave, Tucson, AZ 85721-0065; 
tj arnold@as . arizona.edu 



QSOs, triggered and fueled by galaxy mergers or at least 
interactions, even though there is no evidence to support 
the cla im that mergers are the trigger fo r low- luminosit y 
AGN dFuentes- Williams fc Stocklfl988l: ISchrnlttl 12051 . 
If gas-rich mergers or interactions are the primary trigger 
for these lower-luminosity AGN, there should be higher 
AGN fractions in environments where galaxies have 
an abundant supply of gas and frequent interactions. 
While the cluster environment has very high number 
densities, the galaxies in the centers of rich clusters 
have less cold gas than those in less dense environments 
(e.g. iGiovanelli fc Havnesl [r985). In addition, the high 
pairwise velocity dispersions of cluster members may be 
too large to allow th e formation of boun d pairs inside 
the virial radius (e.g. iGhigna et al"1ll998D . In contrast, 
galaxies in the field have abundant supplies of cold 
gas, but the relatively low galaxy number density may 
counterbalance the abundance of fuel. Between these 
two extremes, the intermediate group environment could 
provide the ideal circumstances for the triggering and 
fueling of low-luminosity AGN in the nearby Universe, 
at least if mergers and interactions make a substantial 
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contribution (although see [MSElisJ L2QQ4j). Galaxies in 
groups have sufficiently modest velocity dispersions, 
numerous neighbors, and available cold gas to trigger 
and fuel AGN. 

Many observations have shown that AGN are rarer 
in clusters when selected via em ission-line diagnostics 
(|Gislerlll978t iDressler et al.lll985t) . Recent observations 
with the the Sloan Digital Sky Survey have shown this 
holds across a wide range in galaxy num ber density, 
specifically for the mo st lu minous AGN dMiller et alJ 
120031 : iKauffmann et alJ 12004 iPopesso fc Biviand I2006T T 
These results have readily been explained with the ar- 
gument outlined above, namely that insufficient cold gas 
reservoirs and high velocity dispersions discourage the 
triggering and fueling of AGN in the cluster environment. 
However, these same observations have indicated that the 
differences between environments fade for less luminous 
AGN, and additionally samples of lower-luminosity AGN 
may vary substantially depending on the selection tech- 
nique employed. The differences in member morphology 
and physical environment between clusters and the field, 
as well as the relati vely unbiased nature of X-ray ob- 
servations, prompted iMartini et al.l (j2006h to conduct a 
survey for X-ray AGN in nearby clusters. They found a 
factor of five higher AGN fraction than previous studies, 
where the AGN were identified to have X-ray luminosi- 
ties Lx > 10 41 erg s _1 in host galaxies more luminous 
than Mr — —20 mag. This dramatic increase was at- 
tributed to the greater sensitivity of X-ray observations 
to lower-luminosity AGN relative to visible-wavelength 
emiss ion-line diagnostics dBaldwin, Phillips, fe Terlevichl 
|1981| ) because the contrast between X-ray emission from 
AGN and other physical processes (low-mass X-ray bi- 
naries, hot gas, and star formation) is higher than in the 
case of the emission-line diagnostics. 

To date only a small number of studies have attempted 
to extend work on X-ray selected AGN to the lower- 
velocity dispersion group environment. The first of these 
was a stud y of op tically-selected groups at z ~ 0.06 by 
iShen et aTl (|2007h . These authors only identified one 
AGN via X-ray selection out of 140 galaxies in eight 
groups , yet found five base d on emission lines. More re- 
cently, Sivako ff et al.1 (|2008l ) compared two, more massive 
groups and four clusters at similar redshifts and found a 
significantly higher X-ray selected AGN fraction in the 
groups compared to the clusters. One of the main goals 
of the present study is to dramatically improve on the 
small group sample in the iSivakoff et al.l (|2008f ) study, 
as well as to provide a larger sample of more massive 
groups than those in the lShen et alJ (|2007l ) study to bet- 
ter span the range of galaxy d ensity from the fi eld to 
clusters. While the groups in the lShen et all p007h sam- 
ple are representative of those found in redshift surveys, 
they are otherwise a fairly heterogeneous sample. In con- 
trast, the X-ray selection of this sample strongly suggests 
that these are all virialized systems. This local sam- 
ple also provides a valuable benchmark for observations 
of the AGN fraction in groups and clusters at higher 
redshifts fe.g. Ueltema et al.l 2007; Silver man et al.ll2009l : 
IMartini. Sivakoff. fc Mulchaeyl I2009T) . For example, at 
z = 0.5 — 1 iGeorgakakis et al.l (|2008f ) find a similar frac- 
ti on of X-ray selecte d AGN in optically-selected groups 
as lShen et al.l (|2007l ) in the local universe. 

Our other motivation is to investigate the role of galaxy 



morphology on the observed AGN fraction, as well as 
on the selection of AGN. Observations of many galaxies 
have shown that there is a strong correlation between 
the mass of the SMBH at the center of a given galaxy 
and the velocity dispersion or luminosity of the host 
galaxy's spheroid component (iFerrarese fe Merrittl 120001 : 
iGebhardt et~aTI l2000l: iMarconi fc Huntl I2003D . These 
relationships imply the coevolution of SMBHs and 
spheroids, and many authors have suggested that 
AGN may actively impact the evol ution of the host 
galaxy by quenching star formation ([Silk fc Reesl 119981: 
Di Matteo. Springel. fc Hernauistl 120051: iHopkins et al.1 
20051: iSpringel. Pi Matteo. fc Hernquistl I2005D . Since 
AGN are a consequence of matter accreting onto the 
central SMBH, and the relation between SMBH mass 
and bulge properties implies a connection in their evo- 
lution, there may be a relation between the in cidence of 
AGN and galaxy morphology. Observations by I Ho et al.l 
(|1997h with the nearby Palom ar Seyfert Survey, and more 
recent work with SDSS (Kauffma nn et aill2003f ). indicate 
that this is the case. These studies find that Seyferts are 
preferentially found in early-type spirals with a signifi- 
cant bulge. One potential interpretation of this result is 
that because early-type galaxies of a given mass will have 
larger SMBHs compared to late- type galaxies of the same 
mass, if their SMBHs both accrete at the same fixed frac- 
tion of the Eddington rate, then the early-type galaxy is 
more likely to be detected in a luminosity-limited sample. 

This is important for comparing galaxy populations 
between groups, clusters, and the field because galaxy 
morphology is obs erved to be a strong function of lo- 
cal galaxy density (Dressier 1980). Furthermore, many 
of the physical processes that are invoked to explain 
the relation between galaxy morphology and density 
may also have relevance for the available fuel supply 
for AGN. These include mergers and ram- pressure strip- 
ping via interactions with the hot ICM (|Gunn fc Gotti 
Il972t lOuilis. Moore, fc Bower]l2000D. evaporation o f the 
cold ISM bv the host ICM (|Cowie fc Songailalll977D . and 
starv ation of new g as that would otherwise replenish the 
ISM ([Larson et alJH$M iBaiogh et al.ll2000h . There are 
thus correlations between morphology and the incidence 
of AGN, and between morphology and galaxy environ- 
ment. It is therefore necessary to take morphology into 
account in order to determine if any variation in AGN 
fraction with environment is directly due to the environ- 
ment itself, or indirectly due to the change in the mix of 
galaxy types with environment. 

To disentangle the effects of morphology and envi- 
ronment we present new observations of rich groups 
and poor cluste r s sele cted from the NORAS sample of 
iBohringer et al.l (|2000l ). These observations and the X- 
ray AGN classification procedure are described in $2] In 
addition, we use spectroscopic data from SDSS to clas- 
sify AGN based on emission-line diagnostics and com- 
pare these AGN to those selected via X-ray emission. 
We then combine t h ese new data with previous work 
([Martini et al.l 120061 : ISivakoff et all |2008) to study the 
incidence of AGN as a function of environment and mor- 
phology. The morphological classification is described 
in detail in §5\ and the AGN fractions are presented 
in 21 The results, including a statistical analysis, are 
described in The final section contains a summary 
of our main results. Throughout this paper we assume 
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H = 70 km s -1 Mpc" 1 . 

2. OBSERVATIONS AND SAMPLE SELECTION 
2.1. XMM-Newton Observations 

Previous studies of X-ray AGN and environment 
have concentrated on rich cluste rs (iMartini et al.l 
2006) and optically-sel ected groups (|Shen et alJ 120071 : 
Gcorgaka kis et all I2008T ) . X-ray detected groups repre- 
sent the intermediate mass-scale between rich clusters 
and poor groups. To study the X-ray AGN population 
in X-ray detected groups, we observed a sample of seven 
low-redshift X-ray groups with the XMM-Newton tele- 
scope. The groups were selected from the NORAS cat- 
alog, which provides a large, uniform sample of X-ray 
bright g roups and clusters fo und in the ROSAT All-Sky 
burvey (jBohringer et aU EToOO) . Our groups were selected 
based on the following three criteria: 1) X-ray luminosi- 
ties between 3 x 10 42 erg s _1 and 3 x 10 43 erg s _1 in the 
0.1-2.4 keV band. This luminosity range was chosen to be 
comparable to the luminosity range of t he intermediate 
redshift X-ray groups in the sample of iMulchaev et ahl 
(2006), potentially allowing a comparison of similar sys- 
tems from redshift zero to z ~ 0.5; 2) Redshifts between 
0.04 and 0.06. This redshift slice was chosen because it 
allows all of the groups to be studied out to approxi- 
mately the virial radius in a single XMM-Newton point- 
ing; 3) Spectroscopic coverage from the SDSS. SDSS 
spectroscopic coverage provides good membership infor- 
mation for each group. This last criterion also allows 
us to calculate an estimate of the AGN fraction based 
on standard emission-line diagnostics (see § 12.2.21 below) 
and to perform morphological fits of the surface bright- 
ness profiles of these objects from the SDSS imaging data 
(see §[3]). The above selection criteria result in a sample 
of 14 X-ray luminous groups and poor clusters. From 
these 14 groups, we selected seven that span the X-ray 
luminosity range of interest. 

The seven X-ray groups were observed by XMM- 
Newton between May 2007 and December 2007. The 
integration times varied from 12 to 40 ksec. All obser- 
vations were obtained in full-frame mode with the thin 
optical blocking filter. The data were processed with ver- 
sion 7.0 of the XMMSAS software. For EMOS data we 
used only patterns 0-12 and apply the #XMMEA_EM 
flag filtering, while for EPN data we used patterns 0- 
4 and set the flag equal to zero. We eliminated peri- 
ods of high b a ckgro und using the method described in 
Ueltema et al.l (|2006| ). This procedure involves filtering 
the data in several energy bands. We begin by applying 
a cut on the high energy (>10 keV) count rate of 0.35 
counts s~ 4 for EMOS data and 1.0 counts s _1 for EPN 
data. We then applied a 3<r clipping to the source-free 
count rate in three energy bands. For this process, we 
used time bins of length 100 seconds. Time bins with 
rates more than 3<r from the mean are then removed un- 
til the mean is stable. Background flaring was very severe 
for one of the groups (RXCJ1225. 2+3213) and resulted 
in no usable data. We therefore eliminated this group 
from our sample. For the remaining six groups, the final 
exposure times were in the range ~ 10 to 23 ksec for the 
EMOS detectors. We note that we only include groups 
in our XMM-Newton analysis for which we can reach a 
lower limit of Lx = 10 41 erg s _1 at a radius of 13'. 



To better constrain the AGN population in groups, 
we supplemented our sample with additional groups 
that have XMM-Newton observations available from the 
archive. These groups were also selected from the NO- 
RAS catalog in the same X-ray luminosity range de- 
scribed above (criterion 1). We also required that these 
groups had spectroscopic coverage in the SDSS (crite- 
rion 3). However, we did not require these groups to be 
in the redshift range 0.04 to 0.06. To insure that a sig- 
nificant fraction of the group members were within the 
XMM-Newton field of view we only considered groups 
that had XMM-Newton coverage out to a radius of at 
least 250 kpc. Using these revised criteria, we added an- 
other five groups to our sample, bringing the total num- 
ber of NORAS-selected groups to eleven. The XMM- 
Newton observations for these additional groups were re- 
duced following the same method described above. The 
final exposure times for all five of these groups is at least 
10 ksec in the EMOS detectors. 

Member ship was determined with the method de- 
scribed in Mulcha ev et al.1 (|2006h and largely based on 
the available SDSS spectroscopy in the fields on these 
groups and clusters. We start with all galaxies located 
within a projected distance of 1 Mpc from the center 
with a recessional velocity within ±3000 km s _1 of the 
mean velocity. We then calculate the velocit y dispersion 
of th e system using the biweight estimator (jBeers et al.1 
1990). Objects with velocities greater than three times 
the velocity dispersion away from the mean are then re- 
moved from the sample and a new mean and dispersion 
are calculated. This process is continued until no further 
objects are removed. Note that although we have esti- 
mated the global properties from galaxies located within 
1 Mpc of the center, our AGN analysis is restricted to 
the smaller radii probed by the XMM-Newton data. 

We supplemented these observations of NORAS groups 
and clusters with other, primarily rich clusters, with red- 
shifts in the range 0.05 < z < 0.08. The det a ils of these 
observations are pre sented in iMartini et al.l (|2006h and 
ISivakoff et al.l (|2008l ). Briefly, these studies are based on 
X-ray observations with Chandra and ground-based im- 
ages and spectroscopy from Las Campanas Observatory 
in Chile. Further details are provided in these two pa- 
pers. Our final list of groups and clusters is provided in 
Table [TJ along with the source of the data for each group 
or cluster. Throughout this study we separate groups 
and clusters based on whether the velocity dispersion is 
greater or less than 500km s _1 . Some of the implica- 
tions of this choice are discussed in §4.11 The division of 
the groups and clusters in Table Q] reflect this threshold 
value. 

2.2. AGN Classification 

While very luminous AGN can be unambiguously iden- 
tified in almost any energy band, AGN become progres- 
sively more challenging to identify at lower luminosi- 
ties when their emission may be equal or even substan- 
tially less than that of their host galaxy. These lower- 
luminosity AGN are important to identify to maximize 
the sample of AGN for demographic studies. It is also 
important to understand the completeness of the AGN 
selection to connect to other studies. At a minimum, 
the completeness should be expressed in terms of lumi- 
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TABLE 1 
Groups and Clusters 



Group/Cluster Name 


Alternate Name 


RA 


Dec 


Redshift 


— Members 


— 




Data Sources 












N 


N* 




[km s —1 ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


A85 




00:41:50.4 


-09:18:11 


0.0554 


109 




53 


993+of 

— OO 


1,3 


A644 




08:17:25.6 


-07:30:45 


0.0701 


75 




40 


952+2o? 

^ — 382 


2,3 


A3128 




03:30:43.8 


-52:31:30 


0.0595 


67 




28 


906 +7 4 


2,3 


RXCJ0110. 0+1358 




01:10:05.5 


+13:58:49 


0.0581 


30 




15 


745+1;! 

' — 64 


1,4 


RXCJ0746. 6+3100 


ZwC10743. 5+3110 


07:46:37.3 


+31:00:49 


0.0579 


23 




16 


719+91 

' — 59 


1,4 


RXCJ1022. 0+3830 




10:22:04.7 


+38:30:43 


0.0544 


36 




18 


710+" 

' xu — 54 


1,4 


A3125 




03:27:17.9 


-53:29:37 


0.0616 


20 




15 


475+ 94 


2,3 


A89B 




00:42:54.6 


-09:13:50 


0.0770 


22 




12 


474+ 155 


1,3 


RXCJ0844.9+4258 




08:44:56.7 


+42:58:54 


0.0550 


13 




9 


343+34 


1,4 


RXCJ1002.6+3241 


ZwC10959.6+3257 


10:02:38.6 


+32:41:58 


0.0505 


33 




9 


454±| 


1,4 


RXCJ1122. 2+6712 




11:22:14.5 


+67:12:46 


0.0553 


22 




8 


223+23 


1,4 


RXCJ1204.4+0154 


MKW4 


12:04:25.6 


+01:54:04 


0.0203 


12 




7 


4951S 


1,4 


RXCJ1223. 1+1037 


NGC4325 


12:23:06.5 


+10:37:26 


0.0255 


4 




2 


334+32 


1,4 


RXCJ1324. 1+1358 


NGC5129 


13:24:11.9 


+13:58:45 


0.0233 


6 




3 


303+| 


1,4 


RXCJ1440.6+0328 


MKW8 


14:40:38.2 


+03:28:25 


0.0269 


15 




9 


449t« 


1,4 


RXCJ1604.9+2355 


AWM4 


16:04:57.0 


+23:55:14 


0.0321 


9 




3 


423+J* 


1,4 


Note. — Properties of all groups and clusters employed 


in this study. 


Columns . 


are: (1) 


Name in the 


original NORAS catalog or 



and 7) number of members more luminous than Mr — — 20 and M% + 1, respectively; (8 ) velocity dispersion; (9) refer ences for data- 
Refer ences are 1: SDSS DR6 images a nd spectroscopy l|A delrrian-McCarthy et al. 2008); 2: 2.5m du Pont Telescope lIMartini et al.l 
2006); 3: Chandra X-ray Observatory jMartini et al.11200 6: Sivak off et al.N2008ri ~4: XMM-Newton, this work. 



nosity in some band, although results are more readily 
compared with theory if they can be expressed in terms 
of bolometric luminosity or accretion rate relative to the 
mass of the black hole. Here we identify AGN via their X- 
ray luminosity, which is estimated to represent on order 
10% of the bolometric luminosity with small scatter (e.g. 
lElvis et"al"1ll994t IMarconi et al.|[2"004f ) and consequently 
it is a reasonable proxy. X-rays also have the advan- 
tage that they are relatively less sensitive to the effects 
of extinction. For the low-luminosity AGN we consider 
here, other physical processes can also produce compara- 
ble X-ray emission from galaxies, so in the first subsec- 
tion below we describe our AGN classification technique 
in detail. The main alternative method to identify AGN 
is via visible- wavelength emission-line ratios and we com- 
pare our X-ray classification to this other method in the 
following subsection (§ I2.2.2|) . 

In addition to careful selection of AGN via either 
method, characterization of how the AGN population 
varies across different environments can be reasonably 
performed with a measurement of the fraction of all 
galaxies of a given morphology that h o st AG N. In pre- 
vious studies (e.g. iMartini et all 120021 l200l the AGN 
fraction was defined as the fraction of galaxies with ab- 
solute magnitude Mr < —20 mag (Vega) that host 
AGN with a broad-band [0.3-8keV] X-ray luminosity of 
Lx > 10 41 erg s _1 . To identify a comparable host lumi- 
nosity range with the NORAS sample, we converted the 
SDSS extinction-corrected, r— band magnitudes (on the 
AB system) to Bessel R— band (on the Vega system) and 
applied a mean k— correction for each group. For both 
of these steps we emp l oyed the software tools described 
by iBlanton fc Roweid ((2007) . We also adopt the evolv- 
ing absolute magnitude threshold of Mr < M R {z) + 1 



introduced by IMartini. Sivakoff. fe Mulchaevl i|2009l ) to 
compare samples across a wide range in redshift, where 
M* R (z) = M R (0) - z, and M* R {z = 0) = -21.92. For the 
present sample the evolution term (z) is negligible and 
the main result is a second threshold approximately one 
magnitude more luminous than the previous magnitude 
cut. This higher threshold is useful because the AGN 
fracti on increases when a higher luminosity threshold is 
used (ISivakoff et al.ll200l . 

2.2.1. X-ray Classification 

We created images in the 0.5 to 8 keV band for each de- 
tector using the flare-cleaned event files spatially binned 
to give 2" pixels. The images for the three detectors 
were combined to form a final image using the SAS task 
EMOSAIC. To identify X-ray sources, we ran the task 
ewavelet on the merged image with a detection thresh- 
old of 5a. The X-ray detections were then compared to 
the known members to determine matches. We restrict 
this analysis to objects within 13' of the field center. 
All sources within 2" of the center of a known mem- 
ber are considered matches; we are motivated to use this 
search rad ius by the la positi onal uncertainty of XMM 
of 1 - 2" (| Watson et al.ll2009f K As is typical of nearby 
X-ray groups and poor clusters, the diffuse X-ray emis- 
sion is centered on or near the brightest galaxy in most 
cases (jMulchaev fc Zabludofl 119981 : [Osmond fc Ponmanl 
l2004h . This makes searching for an AGN component 
difficult for the central galaxies and we have therefore 
excluded these objects from our analysis. This is also 
the case with the analysis of the richer clusters in our 
sample (jMartini et al.ll2006fh 

For each member detected by XMM, we extract a sur- 
face brightness profile to determine the extent of the X- 
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ray emission. We extracted source spectra in circular 
regions extending to where the surface brightness profile 
reaches the background level. Local background spectra 
were extracted from annular regions immediately sur- 
rounding the source. Using a local background of this 
type includes any additional background from the dif- 
fuse intragroup medium at the location of the source. 
Response files (RMFs and ARFs) were constructed for 
the location of the source using the SAS tasks rmfgen 
and ARFGEN, respectively. Source spectra were binned 
to have 25 counts bin -1 . 

All spectral fitting was performed using XSPEC (ver- 
sion 12.3). As noted briefly above, the four main phys- 
ical processes that can produce substantial (Lx > 10 40 
erg s -1 ) broad-band (0.3-8keV) X-ray emission from 
galaxies are AGN, a population of low-mass X-ray bi- 
naries (LMXBs), thermal emission from hot gas, and 
emission from the high-mass X-ray binaries (HMXBs) 
and supernova remnants associated with substantial re- 
cent star formation. We classify an X-ray source as an 
AGN if the observed X-ray luminosity exceeds that ex- 
pected from these other physical processes that can pro- 
duce X-ray emission . Our basic procedure is described in 
ISivakoff et all (l200l. which improves on the earlier pro- 
cedure employed bv lMartini et al.l ()2006l ) . Briefly, we use 
relations between X-ray luminosity, K— band luminosity, 
and star formation rate to determine the expected con- 
tribution from LMXBs (cf. iKim fc Fabbiancj | 2003) . star 
formation (cf. iGrimm. Gilfanov, fc Sunvaevll2003h . and 
halos of hot gas (cf. ISun et al.ll2007f ) and classify a galaxy 
as an AGN if the X-ray luminosity exceeds the expected 
contribution from these other sources of emission. 

Because these XMM observations often have sufficient 
counts for spectral fits, we fit two spectral models to 
better classify AGN when the data are sufficient. These 
models are a single power-law component to represent 
the combined emission of the LMXBs and any AGN 
component and a thermal component to represent any 
emission from hot gas. We then estimate the X-ray bi- 
nary emission expected from the if— band luminosity of 
the galaxy and consider any excess emission from the 
power-law component to be due to an AGN. In all our 
fits the neutral hydroge n column densi t y is fi xed at the 
Galactic value given in Kal berla et all ([2005h . We also 
fix the power-law index (r) to 1.7. For the thermal com- 
ponent, we use the MEKAL model in XSPEC with the 
abundance fixed at 0.8 solar. We simultaneously fit the 
spectra from all three EPIC detectors. For galaxies with 
at least several hundred counts, it is usually possible to 
constrain both the thermal and powerlaw components. 
In some cases, only one component is required to pro- 
duce an adequate fit (i.e. the normalization of the sec- 
ond component is consistent with zero). For galaxies 
with a small number of counts, it is not possible to dis- 
tinguish between the possible spectral models. For these 
objects, we have estimated the X-ray luminosity assum- 
ing a power-law model alone. The resulting luminosities 
of the power-law and thermal components for each galaxy 
are given in Table El We note that if a thermal model is 
assumed for the cases where the spectral model cannot 
be determined, the resulting luminosities would be lower 
by a factor of approximately two. 

The additional spectral information available for many 
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Fig. 1. — The broad-band X-ray luminosity Lx,B vs. the near-IR 
luminosity, Lk s , tot for the galaxies displayed in Table [2] The top 
panel shows the broad-band X-ray luminosity Lx b 01 the best-fit 
power-law component for all of these eight galaxies compared to 
the K— band luminosity. Here we compare to the expe cted relation 
for LMXBs {dotted line) from IKim fc F abbiano (2004) and classify 
four of eight galaxies as X-ray AGN (stars). The bottom panel 
shows the relation for sources with insufficient counts for spectral 
modeling. Here we assume a T = 1.7 power law and compare 
with the sum (solid lin e) of the LMXB relation and a thermal 
model (dashed line) from Sun et al. ( 2007). Five of the sources are 
classified as AGN, although one is below 10 41 erg s — 1 and is not 
included in the statistical analysis. In both panels the thicker line 
is the relation and the thinner lines are ±lcr uncertainties. Other 
galaxies are classified as inactive (filled circles, see Table [2l . 

of these sources better constrains the nature of the X- 
ray emission. Specifically, eight of the 14 X-ray sources 
associated with members have sufficient signal-to-noise 
ratio to determine if the X-ray emission was best-fit by a 
power-law model, thermal m odel, or both. This im proves 
the classification over, e.g. ISivakoff et all |2008) as we 
can then compare the luminosity of the best-fit power- 
law model to just the expected emission from LMXBs, 
HMXBs, and an AGN component, and exclude the ther- 
mal model because of its different spectrum. The top 
panel of Figure [1] shows the broad-band X-ray luminos- 
ity of the best-fit power-law component for all of these 
eight galaxies c ompared to the K— b and luminosity. Note 
that whereas in lSivakoff et all (|2008T ) we used the 2MASS 
K20 magnitude, here we employ the Ktotai magnitude as 
this is a better match to the aperture used for the X-ray 
photometry. The figure also shows the expected rela- 
tionship between X-ray luminosity and K— ban d lumi- 
nosity for LMXBs from IKim fc Fabbianol (|2004[ ). where 
the thicker line is the relation and the thinner lines are 
±1(T uncertainties. Four of the eight X-ray sources fall 
on the LMXB relation and we classify these galaxies as 
inactive (see Tabled]). The remaining four are at least 
2cr more X-ray luminous than would be expected from 
LMXBs alone and we therefore classify these galaxies 
as X-ray AGN. We note that these classifications are 
the same as we would have assigned based on our pre- 
yious approach with an Lx = 10 41 er g s -1 threshold 
(|Martini et al.ll2(M ISivakoff et al.ll200l) . 

The other six galaxies had sufficiently faint X-ray emis- 
sion that we were unable to accurately model their X-ray 
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TABLE 2 
X-ray Detected Group Members 



Galaxy 




Lk 


Model 


L Xf 

10 4({ 

(0.3- 


>owcrlaw 

erg s _1 
8.0 keV) 


Lxthermal 

10 40 erg s" 1 
(0.3-8.0 keV) 


Class 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 



2MASXJ07463295+3101213 


-22 


.6 


11 


.12 


Po 1 




18.6+JV 




AGN 


2MASXJ07462331+3101183 


-22 


.2 


11 


.36 


Po 1 








inactive 


2MASXJ08445063+4302479* 


-23 


.1 


11 


.67 


Po 




24.iit-;° 




AGN 


2MASXJ10230356+3838176 


-21 


.0 


10 


.70 


Po 1 




79 9 +31.4 
,z - z -30.8 




AGN 


2MASXJ10223745+3834447 


-23 


.2 


11 


.69 


Po+ 


Th 


104.9^-8 


33.0lf 5 5 


AGN 


2MASXJ10220069+3829145 


-21 


.9 


11 


.03 


Po 1 




OO Q + 8.8 

§ 2 




AGN 


2MASXJ11231618+6706308 


-22 


.1 


11 


.22 


Po 1 




le.e+^g 7 




AGN 


2MASXJ11221610+6711219 


-21 


.5 


10 


.96 


Po 




i4.ol^; 3 5 




AGN 


2MASXJ11223691+6710171 n 


-21 


.6 


11 


.03 


Po 1 




8.5l« 2 9 




AGN 


SDSSJ112333.56+671109.9 


-20 


.1 






Po+ 


Th 




14.9+fs 1 


AGN 


2MASXJ12043806+0147156 


-22 


.5 


11 


.33 


Po+ 


Th 


7 n+3.7 
'• z -2.0 




inactive 


2MASXJ12225772+1032540 


-21 


.3 


11 


.08 


Po+ 


Th 






inactive 


2MASXJ13242889+1405332 


-22 


.3 


11 


.43 


Po 








inactive 


2MASXJ14403793+0322375 


-22 


.3 


11 


.38 


Po 




"•-■■-l.g 




inactive 



Note. — X-ray measurements and classifications. Columns arc: (1) Galaxy name; (2) Host galaxy 
-R-band absolute magnitude; (3) Total if-band luminosity from 2MASS (for SDSSJ112333. 56+671109. 9 
K — 14 mag was assumed based on the colors of other group members); (4) Model fit to the X-ray 
data where Po is a power-law fit, Th is a thermal model, and Po 1 indicates a power-law was assumed (see 
o! !2.2.H ; (5) X-ray luminosity of the power law component; (6) X-ray luminosity of the thermal component; 
(7) classification of the galaxy as either an AGN or as inactive. Lx,po and Lx.th arc in units of 10 40 
erg s — 1 and are broad-band (0.3-8kcV) measurements. The * superscript in Column 1 denotes the single 
galaxy in the NORAS sample that is classified as an AGN based both on its X-ray properties and its 
emission lines. The n superscript refers to an AGN below our Lx — 10 41 erg s" 1 limit and thus not 
included in the sample statistics. 



spectra. In these cas es we followed the procedure used by 
ISivakoff et al.1 (|2008l ) and assumed a T = 1.7 model and 
measured the X-ray luminosity of that model. These six 
sources are shown in the bottom panel of Figure [T] along 
with the same LMXB relation shown in the top panel, a 
relation between Lx a nd L k for t herma l emission from 
hot gas adapted from iSun et all (|2007D . and the sum 
of these two relations, where again the inactive galaxy 
relation is represented by the thicker line and the ±lcr 
uncertainties are represented by thinner lines . The hot 
gas re lation is modified from that presented bv lSun et al.l 
(|2007f ) because their measurements were in the soft band 
(0.5-2keV) and ours are broad-band (0.3-8keV) measure- 
ments. We therefore multiplied their X-ray luminosity by 
the flux ratio of a kT = 0.7keV thermal bremsstrahlung 
model in the broad band and soft band. In practice 
none of the sources in the lower panel are sufficiently 
luminous in the K— band that we expect a substantial 
thermal component. Five of the six sources are above 
the relation by at least 2a and we classify these galaxies 
as X-ray AGN, although one of these is not included 
in our statistical analysis because it has Lx < 10 41 
erg s _1 . We also note that none of these 14 galaxies 
appears to have sufficient star formation to contribute 
significantly to the X-ray luminosity. To check this we 
estimated sta r formation rates for each galaxy from the 
Ha flux (e.g. lKennicurt1ll998t iBrinchmann et al.ll2004h 
and then calculated the expected X -ray luminosity from 
IGrimm, Gilfanov. fc Sunvaevl ((2003) . In all cases the es- 
timated X-ray luminosity due to star formation was at 
least an order of magnitude below the observed luminos- 
ity. For the five groups and clusters in common with 



ISivakoff et all (12001 we retain the same AGN classifica- 
tions presented in that paper. 

2.2.2. Emission-Line Classification 

The traditional technique to identify low-luminosity 
AGN in emission-line galaxies is through use of line ra- 
tios such as [OIII1A500 7/H/3 and [NII1A6584/Ha on a 
"BPT Diagram" (jBaldwin. Phillips, fc Terlevich] [l98lh . 
We have measurements of these four emission lines for 
most (349) of the galaxies in the sample from the NORAS 
catalog from the MPA-JHU galaxy catalog for the SDSS 
Data Release 7 2 . We use the emission-line measurements 
and errors calculated by the MPA-JHU group, which 
have been corrected for stellar absorption, and then cal- 
culate emission line ratios. If the line flux measurements 
for [NII]A6584/Ha are larger than 3 times the error as- 
sociated with the measurements, we keep the data for 
analysis. We are less conservative for the [OIII] A5007/H/? 
ratio because galaxies with precise [NII]A6584/Ho; mea- 
surements can be unambiguously identified as AGN even 
if the [OIII1A5007 /H/3 measurements are less certain (e.g. 
iShen et aTTl2007t ). In addition, edge-on galaxies might 
obscure the bluer [OIII] and H/3 lines, leading to a mis- 
classification. We hope to keep these possibly obscured 
but still unambiguous AGN in our sample by being less 
stringent with our condition on the [OIII]A5007/H/3 mea- 
surement errors. 

Of the 349 galaxies in the MPA-JHU galaxy catalog, 
116 have sufficiently bright spectral lines to place on the 
BPT Diagram shown in Figure [2j This figure illustrates 

2 http:/ /www. mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 2.— A BPT Diagram (Baldwin, Phillips, fc TerlevichllT98lh 
for all galaxies with sufficiently bright emission-lines in our NORAS 
sample. Inactive galaxies (open circles) are separated based on 
whether Mr > —20 mag (small open circles), —20 > Mr > M R +1 
(medium open circles), or have Mr < M R + 1 (large open circles). 
The X-ray AGN (filled circles) and BPT AGN (stars) are similarly 
scaled based on host galaxy luminosity and all AGN have error 
bars. Galaxies are classified as BPT AGN if they fall above the 
I Kewlev et ail 120011) cri terion (dot- dashed line). The criterion of 
Kauffmann ct al. (2003) is also shown (dashed line). 

the values of [NII]A6584/Ha and [OIII]A5007/H/3 for 
these 116 galaxies. We identify galaxies as "BPT AGN" 
if these two line ratios place the m above the threshold 
suggested by IKewlev et al.l (|2001[ ) as the maximum limit 
that could be produced by extreme starburst activity. 
This is a very conservative limit and many confirmed 
AGN do not meet this criterion. For reference, we also 
show the less conservative threshold of iKauffmann et "aTl 
(2003). Of the 116 galaxies with sufficient emission- line 
detections to fall on this diagram, 14 were identified as 
BPT AGN. A total of 199 out of the 349 galaxies with 
SDSS spectra also have M R < -20 mag, so the BPT 
AGN fraction for this sample is 14/199 = 0.070±^ig (all 
BPT AGN have Mr < — 20 mag). For comparison, six of 
these 199 galaxies or ~ 3% were identified as X-ray AGN 
(note that not all X-ray AGN had sufficient spectral-line 
data). It is also striking that only one galaxy is classified 
as both an X-ray AGN and a BPT AGN. The samples are 
nearly disjoint. Further information about the galaxies 
identified as BPT AGN is provided in Table [3] 

3. MORPHOLOGICAL ANALYSIS 

Though galaxy morphology has been used extensively 
to study galaxy evolution, until recently this property 
was determined by eye. This process, though largely 
repeatable, lacks quantitative robustness and is a pro- 
tracted process for large numbers of objects. The 
alternative — measuring morphology in an automated 
fashion — is not simple to implement. Only in the last 
decade, with the development of large- format, linear de- 
tectors and substantial computational resources, has it 
become commonplace to classify galaxies using quanti- 
tative and repeatable techniques. Datasets consisting of 
large numbers of galaxies have made visual inspection 
intractable as a method for determining morphology (al- 



though there have been novel approaches f or morpholog- 
ical i dentification via visual inspection; see lLintott et al.l 
2008), while abundant computational resources have 
made large-scale, quantitative analyses more feasible. 

Various methods to ascertain morphol- 
ogy quantitatively exist in t he literature 

(Conselice. Bcrsha dv. fc Jangrenl [20001: ISimardl [2002; 
iGoto et alj 120031 : lLotz et al.ll2004[ ) , and we chose to use 
the galaxy fitting code GALFIT (|Peng et all I2002D to 
measure the morphological properties for galaxies in 
our sample. One of our main m otivations for th e choic e 
of GALFIT was the work of lHaeussler et all (120071) . 
who compared GALFIT to GIM2D (jSimardJ |2002j) and 
concluded that GALFIT has advantages in its ability 
to simultaneously fit neighboring galaxies in a crowded 
field and benefits from a dramatic increase in execution 
speed. GALFIT is designed to extract structural 
components from galaxies by fitting two-dimensional 
light profiles with an arbitrary number of parametric 
functions that are suitable for describing the surface 
brightness distribution of galaxies. Although the code 
was authored to fit subtle structures of well-resolved 
galaxies with many-component models simultaneously, 
it is also effective at handling large numbers of galaxies 
imaged at lower resolution by fitting their surface 
brightness profiles with relatively simple models. We 
utilize the latter capability in our analysis. GALFIT 
takes as input a simple text file and is very customizable, 
allowing easy extension via a wrapper script. A final 
advantage is that GALFIT can use a variety of analytic 
functions singly or simultaneously, including the Sersic 
profile. In the next subsections we describe the fitting 
procedure in more detail, including how we parametrize 
galaxy morphology, and present the results of our fits of 
the X-ray and BPT AGN. We fit models to all of the 
galaxies in our sample that have imaging data. 

3.1. Model Fits 

We have adopted two parameters to quantify galaxy 
morphology: the value of the Sersic index n and the 
bulge-to- total flux ratio B/T. The Sersic profile is de- 
fined to be: 



S(r) = E e exp 




- 1 



(1) 



where E e is the surface brightness (flux per unit area) 
at the effective radius r e , r e is the half-light radius, and 
the eponymous index n characterizes the shape of the 
light profile. The parameter n is set by the constraint 
that r e is the half-light radius. The Sersic index n has 
been commonly used to separate early-type and late-type 
galaxies in the literature. Studies based on SDSS imaging 
data find that n = 2.5 is a reasonable po i nt to distinguish 
these two types (e.g. iBell et al.l 12004a!: iMcIntosh et alj 
l2005f ) and we too adopt n > 2.5 to identify early- type 
galaxies. 

The bulge-to-total flux ratio is measured from a clas- 
sic bulge-disk decomposition. Here we fit a bulge compo- 
nent with a ide Vaucouleursl (|1948l ) r 1 / 4 profile and a disk 
component with an exponential surface brightness pro- 
file. The fraction of the total flux in the bulge component 
relative to the total (bulge + disk) flux then provides the 
ratio B/T. Throughout this work we will primarily rely 
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TABLE 3 
BPT AGN 



Galaxy 


M R 


log([OIII]A5007/H/3) 


log([NII]A6584/Ha) 


(1) 


(2) 


(3) 


(4) 


oi\/T A qv Tm nofinnoxi q^qi 

ZWli\oA.JUlUyDyUZ-r loOolDO 


on q 


a 1 qo 4- n non 


U.14U It U.1DD 


cnGQTmnnc;7 oq i i/inQon i 
olJooJUlUyo ( .oo-+-14UozU. 1 


on a 


-u.uoz ± u.uyy 


U.ooz it u.zyz 


SDSSJ011021. 57+135421.4 


-20.2 


0.180 ± 0.119 


0.501 ± 0.230 


2MASXJ07462331+3101183 


-22.2 


0.164 ± 0.049 


0.285 ± 0.082 


2MASXJ07470054+3058205 


-21.7 


0.129 ± 0.050 


0.207 ± 0.170 


2MASXJ08445063+4302479 * 


-23.1 


0.032 ± 0.084 


0.061 ± 0.078 


SDSSJ100311. 10+323511.3 


-20.2 


0.304 ± 0.062 


0.362 ± 0.361 


2MASXJ10213991+3831195 


-21.5 


-0.038 ± 0.119 


0.278 ± 0.128 


2MASXJ11223691+6710171 (VIIZw392) 


-22.4 


0.168 ± 0.020 


0.005 ± 0.046 


SDSSJ112425. 38+671940.0 


-20.3 


0.115 ± 0.045 


0.281 ± 0.092 


2MASXJ12041899+015054 (CGCG013-058) 


-21.2 


0.129 ± 0.072 


0.064 ± 0.146 


2MASXJ12230667+1037170 (NGC4325) 


-22.5 


0.151 ± 0.027 


-0.251 ± 0.171 


2MASXJ12225772+1032540 


-21.3 


-0.082 ± 0.152 


0.343 ± 0.106 


2MASXJ13241000+1358351 (NGC5129) 


-23.1 


0.162 ± 0.084 


0.086 ± 0.120 



Note. — Galaxies identified as AGN on the BPT Diagram (Baldwin, Phillips, <fe Tcrlcvich 1981) shown 
in Figure [2] Columns arc: (1) The name (alternate name) of the galaxy; (2) A-/r; (3) calculated value 
of log [OIIIJA5007/H/3 from the MPA-JHU database; (4) calculated value of log [Nil] A6584/Hq from the 
MPA-JHU database. The * superscript in Column 1 denotes the only galaxy that is classified as an AGN 
based both on its X-ray properties and its visible-wavelength emission lines. 



on the Sersic index to classify galaxies, but we will use 
the B/T as a consistency check on our results. Note also 
that n = 4 is equivalent to the r 1 / 4 profile and n = 1 is 
equivalent to an exponential disk profile. 

When available we also compared our calculated mor- 
phological parameters to the SDSS f racDeV parameter, 
which is calculated by the SDSS pipeline and serves as 
another quantitative measure of galaxy morphology. The 
f racDeV parameter is very similar to our bulge-to-total 
flux decomposition. It is obtained by fitting the sur- 
face brightness profile of a galaxy with exponential and 
de Vaucouleurs components and then keeping the frac - 
tional contribution of the latter. iBernardi et alj (2005) 
identified f racDeV > 0.8 to identify early-type galaxies. 

As noted above, GALFIT accepts an input text file 
that specifies the initial conditions for models and other 
options and parameters related to the fit. GALFIT can 
accommodate as many independent models as the user 
desires, limited only by computational resources. It is 
also possible to simultaneously fit adjacent or blended 
objects in addition to the target to remove potential 
contamination and obtain a more robust fit. GALFIT 
convolves the model with a point spread function (PSF) 
supplied by the user, subtracts the convolved model from 
the input image, and computes the reduced chi squared, 



ni ny 

EE 



(flux^ - model^y) 



N dol 



(2) 



x—l y—1 



'x,y 



where Adof is the number of degrees of freedom in 
the model, flux and model are the pixel values of the 
original image and analytic model, respectively, and 
a\ is the error in each pixel. GALFIT minimizes 

\t using a Levenberg-Marquardt algorithm, a downhill- 
gradient type algorithm suited for searching large pa- 
rameter spaces quickly. Additional and optional input 
includes a bad pixel map specifying which pixels should 
be excluded from the xt calculation (i.e. masked out) 
and initial guesses for the many free parameters, includ- 



ing the astrometric and morphological quantities of the 
target. 

GALFIT is useful in its extensibility and we took ad- 
vantage of this by creating a wrapper script in Python 
and an algorithm to fit many target objects with little 
to no user interaction. The results are comparable to a 
human user fine tuning GALFIT input parameters until 
an ideal fit is obtained. The input to this process is a list 
of astrometric coordinates of targets and the FITS im- 
ages in which these targets are imaged. For each target 
the script determines if it lies within the boundaries of 
a given image and obtains initial morphological param- 
eters (e.g. a measure of the galaxy's radius, magnitude, 
ellipticity, object position, and posi tion angle) with SEx- 
TRACTOR (|Bertin fc Arn outs 1996|). At this stage we re- 
tain information about all detected objects within some 
arbitrary number (found using trial and error) of effec- 
tive radii from the target. The fitting region supplied 
to GALFIT is determined in a similar way. Based on 
the parameters of the objects in the field of view, we ei- 
ther add the pixels associated with the object to the bad 
pixel file, masking them out and removing them from the 
xt calculation and fitting procedure, or fit the object in 
addition to the target. This discriminatory algorithm 
compares object brightness to the target brightness (e.g. 
an object very much dimmer than the target will likely 
be masked out rather than fit) and the distance from the 
target (e.g. an object separated by many target galaxy- 
radii will likely be masked out rather than fit). In this 
way we only fit additional objects if they are likely to 
contaminate the fit of our target, and thus its morpho- 
logical properties. Fitting a superfluous number of ob- 
jects is computationally wasteful and complicates finding 
a unique minimum in xt space. 

Our wrapper script generates a GALFIT input file for 
both a single component Sersic profile fit and a two com- 
ponent Sersic profile + Exponential Disk fit. Several it- 
erations of this dual method fitting occur if the resultant 
xf, values for the two parameterizations differ by more 
than a small amount, with the previous fit parameters 
used for subsequent iterations, to insure that the solu- 
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tion obtained is not the result of the minimization algo- 
rithm getting lost in a local minimum. This propensity 
to get lost in a local minimum and the question of the 
uniqueness of a multi-component solution in large pa- 
ram eters spaces is an issue for algorithms like GALFIT, 
and iPeng et al.l (|2002f ) addresses this question in depth 
(see their § 3.3). 

Another requirement for GALFIT is an estimate of the 
PSF shape. These were created manually for each im- 
age by approximating them with a Gaussian generated 
with the standard gauss task in IRAF. Effective radii of 
the PSFs were estimated by examining several stars in 
each image. We determined that small variations in the 
radius of the PSF negligibly affected the parameters of 
the model, confirming the robustness of this approximate 
method of PSF creation. 

When completed, the wrapper saves: 1) a postage 
stamp sized image of the original galaxy (i.e. a cropped 
section of the wide-field FITS image where the target 
fills the frame); 2) an image of the same size that dis- 
plays the model generated by GALFIT convolved with 
the user-defined PSF; 3) a residual image: the original 
data less the convolved model. This third image should, 
if the model is perfect, show a scene identical to the input 
image, with the exception that the region of sky previ- 
ously occupied by the target galaxy should be indistin- 
guishable from noise (although somewhat increased noise 
due to the contribution of the subtracted object). Unsur- 
prisingly, this ideal case is rare (though not nonexistent). 
These three images are saved for each fit technique (sin- 
gle and two-fit methods) for quick visual inspection. We 
also retain the fit log files generated by GALFIT, FITS 
images that contain the GALFIT models generated by 
the algorithm, and a text file with the relevant morpho- 
logical parameters and object information for later use. 
The script that automates these tasks works for an arbi- 
trarily large set of input targets. 

3.2. Fit Results 

The three-image results for all of the X-ray AGN are 
shown in Figure while the results for the BPT AGN 
are shown in Figure 3J The model parameters associated 
with all of these fits are listed in Tables Q] and [5l Only one 
object falls in both samples and consequently appears in 
both figures and tables. While these figures demonstrate 
our results on AGN, described further below in H4.2[ they 
also are fairly representative of our morphological fits to 
the inactive group and cluster galaxies. 

One common challenge to most fits is crowding. Fig- 
ure [31 Panel H shows a galaxy in the cluster A644 where 
many neighbors have been removed by subtraction and 
only the target was fit with GALFIT. Our algorithm de- 
termined that all neighbors were sufficiently far or faint 
enough that they did not interfere with the fitting proce- 
dure. To obtain a low value of xt> the pixels from these 
other objects in the field of view were simply masked 
out. Figure [3l Panel D depicts a similar but slightly 
more challenging case. Here the target is blended with 
a bright, nearby object. In this case multiple objects 
are fit to obtain a reasonable model for the target. As 
before, objects yet further away are simply masked out. 
Panels A & B present additional examples of this case. 
A final case is illustrated by Figure [H Panel M. In this 
case the model is simply inadequate because the galaxy's 



morphology is more complicated than our simple mod- 
els. This galaxy appears to be both blended with other 
objects and morphologically disturbed. 

The algorithm used by GALFIT is optimized for speed 
and is searching a very large parameter space; conse- 
quently the probability of getting lost in a local minimum 
is non-negligible. Though it is difficult to be certain that 
this has not happened for most cases, we employ sev- 
eral techniques to guard against this possibility. First, 
we run GALFIT on each target several times and up- 
date the initial fit parameters if the resultant xt differs 
by more than several percent between two runs. Second, 
we fit two distinct sets of models to each galaxy in sepa- 
rate runs: the Sersic index in a single component fit and 
the bulge-to-total flux ratio in a two component model 
fit. These quantities are correlated for most galaxies (see 
Figured]) and we carefully reexamine egregious outliers. 
Finally, we save output images, as shown in Figures [3] 
and[U and inspect these to identify poor fits. 

Stubborn objects that are not well fit by our procedure 
persist, though they are relatively fe w. As noted above , 
GALFIT was originally developed bv lPeng et all (|20Q2f ) 
with the purpose of decomposing the complex structure 
of well resolved galaxies. We instead use it to do simple 
one- or, at most, two-model fits. For our typical resolu- 
tion and galaxy size this is not a problem. However, some 
of the galaxies in our images are so well resolved that a 
simple Sersic profile or de Vaucouleurs bulge plus expo- 
nential disk fit is insufficient to adequately describe the 
morphology. Structures such as bars or prominent spiral 
arms are sometimes fit rather than the more averaged 
profile of the galaxy that may have resulted from a less 
well-resolved image of the same object. Also, we occa- 
sionally find objects that are blended or coincident with 
the target and cannot be simultaneously with the target, 
such as the example of Panel M of Figure [4] mentioned 
above. In general, blending is a particularly challenging 
problem if SExtractor fails to find the blended object 
as a distinct source. 

In addition to comparison of the correlation between 
the Sersic index and the bulge-to-total flux ratio shown in 
Figure El we also compare to the SDSS f racDeV parame- 
ter when available. We find that these two parameters do 
not correlate well, while our Sersic index measurements 
do correlate with the fracDeV parameter, as shown in 
Figure [U although the scatter is significant. One source 
of scatter is that the f racDeV parameter saturates while 
we still measure a substantial range of Sersic index. We 
fit a line to this relation, although excluded points with 
fracDeV = 1, and found that fracDeV = 0.66 corre- 
sponds to n = 2.5. This is reasonably consistent with 
the value of f racDeV = 0.8 used in the literature to iden- 
tify early- type galaxies (e.g.,Bcrnardi et al. 200i). As for 
Figure El we visually inspected all outliers on this figure 
to check the goodness of fit. 

We performed this morphological analysis on all of the 
confirmed me mber galaxies and a dditional groups and 
clusters from (|Sivakoff et all 120081 ). We used these data 
to identify all of the early-type galaxies in each group 
and cluster with Mr < —20 mag and Mr < M R + 1 and 
calculate the early-type galaxy fraction. These results 
are presented in Tabled which includes the total number 
of galaxies and the number successfully fit. Typically 
these numbers agree, except for rare instances when some 
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Fig. 3. — Data, GALFIT mo dels, and residuals of all the X-ray AGN for our entire sample. For details of the AGN identification see 



£12.2.11 Martini ct al. (2006), and Sivakoff et al. (2008). The morphological parameters associated with these fits are listed in Table H The 
star symbol, next to Panel O, denotes the only galaxy in our sample identified as both a BPT AGN and an X-ray AGN. This object is also 
in Panel F in Figure [4] 



galaxies were not fit successfully. 

4. AGN FRACTIONS 

The next step of our analysis is to combine the AGN 
classifications from H2.'2\ and the morphology fits from 
|J3] to measure the AGN fraction as a function of envi- 
ronment and determine if there is any variation with 
host galaxy morphology. This analysis is described in 
the first subsection below. In addition, in §2.21 above we 
demonstrated that the X-ray AGN and BPT AGN were 
nearly disjoint populations. In the following subsection 
we compare the morphologies of these two populations. 
Throughout we calculate AGN fractions for absolute 



magnitude limits of Mr < — 20 mag and Mr < M R + 1 
and in all cases the AGN fraction is defined to be the 
number of AGN divided by the total number galaxies 
above a given absolute magnitude limit. All errorbars 
are derived from Poisson and binomi al statistics an d are 
single-sided, ler confidence intervals (|Gehrelslll986l ). 

4.1. X-ray AGN 

Table [7] provides the number of X-ray identified AGN 
in each group or cluster, the X-ray AGN fraction, and the 
X-ray AGN fraction with early-type host galaxies. These 
results are also illustrated graphically in Figures [7] and 
although groups and clusters with no AGN are not shown 
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Fig. 4. — Data, GALFIT models, and residuals of all the BPT AGN for our entire sample. For details of AGN identification see § 12.2.21 
The morphological parameters associated with these fits are listed in Table [5] The star symbol, next to Panel F, denotes the only galaxy 
in our sample identified as both a BPT AGN and an X-ray AGN. This object is also in Panel O in Figure [3] 



TABLE 4 

Fit Properties from Figure [3] — X-ray AGN 



Panel 


Cluster / Group 


Object Name 


Sersic Index 




fbulge / ftotal 


A U 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


A 


A3125 


J032723.4-532535.5 


3.18 


1.965 


0.930 


1.875 


B 


A3125 


J032725.3-532506.6 


2.27 


1.635 


0.768 


1.638 


C 


A3125 


J032705.1-532140.9 


2.63 


1.739 


0.724 


1.630 


D 


A3128 


J033039.2-523205.7 


1.24 


1.300 


0.109 


1.236 


E 


A3128 


J032941.4-522935.7 


1.36 


1.182 


0.335 


1.178 


F 


A3128 


J033051. 0-523031. 2 


4.53 


1.528 


0.964 


1.652 


G 


A3128 


J033017.3-523408.9 


1.75 


1.585 


0.458 


1.576 


H 


A644 


J081748.1-073731.7 


2.16 


9.746 


0.575 


7.099 


1 


A644 


J081739.5-073309.0 


12.0 


10.58 


0.371 


7.762 


J 


A85 


J004311. 6-093816.1 


3.06 


1.583 


0.400 


1.568 


K 


A85 


J004130.3-091545.9 


3.57 


1.712 


0.085 


1.606 


L 


A89B 


J004314.1-092145.2 


2.89 


1.538 


0.545 


1.532 


M 


RXCJ1122.2+6712 


SDSSJ12333.56+671109.9 


1.56 


1.498 


0.053 


1.495 


N 


RXCJ0746.6+3100 


2MASXJ07463295+3101213 


4.28 


2.044 


1.000 


2.038 


0* 


RXCJ0844.9+4258 


2MASXJ08445063+4302479 (CGCG208-041) 


4.59 


1.551 


0.812 


1.552 


P 


RXCJ1022. 0+3830 


2MASXJ10220069+3829145 


2.37 


6.604 


0.808 


6.556 


Q 


RXCJ1022. 0+3830 


2MASXJ10223745+3834447 (NGC 3219) 


3.49 


6.268 


0.851 


6.268 


R 


RXCJ1022. 0+3830 


2MASXJ10230356+3838176 


9.38 


6.672 


0.534 


6.672 


S 


RXCJ1122. 2+6712 


2MASXJ11221610+6711219 


6.88 


1.625 


0.355 


1.587 




RXCJ1122.2+6712 


2MASXJ11223691+6710171 (VIIZw394) 


3.35 


1.550 


0.894 


1.547 


u 


RXCJ1122. 2+6712 


2MASXJ11231618+6706308 


4.60 


1.520 


1.000 


1.524 


V 


A89B 


J004300.63-091346.4 


5.97 


1.544 


0.999 


1.569 


Note. 


— GALFIT output parameters from fits to all X-ray AGN. The * superscript in Column '. 


identifies the only galaxy in our 



sample that is classified as both a BPT AGN and an X-ray AGN. This object is Object F in Tablc[5] The 71 superscript refers to an 
AGN below our L x = 10 41 erg s -1 limit. 
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TABLE 5 

Fit Properties from Figure |4] — BPT AGN 



Panel 


Cluster / Group 


Object Name 


Scrsic Index 


xl 


fbulgc / ftotal 


xl 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


A 


RXCJ0110.0+1358 


2MASXJ01095902+1358155 


2.99 


1.625 


0.932 


1.648 


B 


RXCJ0110.0+1358 


SDSSJ010957.88+140320.1 


2.20 


1.567 


0.398 


1.564 


C 


RXCJ0110.0+1358 


SDSSJ011021.57+135421.4 


1.96 


1.523 


0.562 


1.524 


D 


RXCJ0746. 6+3100 


2MASXJ07462331+3101183 


2.55 


2.221 


0.509 


2.112 


E 


RXCJ0746. 6+3100 


2MASXJ07470054+3058205 


2.96 


1.604 


0.613 


1.582 


F* 


RXCJ0844.9+4258 


2MASXJ08445063+4302479 (CGCG208-041) 


4.59 


1.551 


0.812 


1.552 


G 


RXCJ1002. 6+3241 


SDSSJ100311. 10+323511.3 


1.51 


1.424 


0.472 


1.427 


H 


RXCJ1022. 0+3830 


2MASXJ10213991+3831195 


1.13 


6.372 


0.022 


6.372 


I 


RXCJ1122. 2+6712 


2MASXJ11221537+671318 (VIIZw392) 


4.86 


1.595 


0.667 


1.582 


.1 


RXCJ1122. 2+6712 


SDSSJ112425.38+671940.0 


1.69 


1.531 


0.509 


1.532 


K 


RXCJ1204.4+0154 


2MASXJ12041899+015054 (CGCG013-058) 


2.17 


1.798 


0.726 


1.747 


L 


RXCJ1223. 1+1037 


2MASXJ12230667+1037170 (NGC4325) 


2.42 


1.627 


0.884 


1.579 


M 


RXCJ1223. 1+1037 


2MASXJ12225772+1032540 


3.96 


2.313 


0.977 


2.316 


N 


RXCJ1324. 1+1358 


2MASXJ13241000+1358351 (NGC5129) 


4.22 


1.805 


0.792 


1.762 



Note. — GALFIT output parameters from fits to all BPT AGN. The * superscript in Column 1 identifies the only galaxy in our 
sample that is classified as both a BPT AGN and an X-ray AGN. This object is Object O in Table [4] 



TABLE 6 
Morphological Results 



Cluster / Group 




N m - 


/n>2.5 








iV fit . - 


/n>2.5 




- N* , 

total 






All 


n > 2.5 




Raw 


Corr 


n > 2.5 All 


n > 2.5 


Raw 


Corr 


n > 2.5 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 


(11) 


(12) 


(13) 



A85 


104 


76 


0. 


.731; 


A644 


15 


7 


0. 


.467; 


A3128 


54 


21 


0. 


389; 


RXCJ0110.0+1358 


30 


19 


0. 


633; 


RXCJ0746.6+3100 


23 


17 


0. 


.739; 


RXCJ1022.0+3830 


36 


11 


0. 


306; 


All Clusters 


262 


151 


0. 


576 



A3125 


IS 


10 


0. 


556; 


A89B 


14 


7 


0. 


500; 


RXCJ0844.9+4258 


13 


10 


0. 


.769; 


RXCJ1002.6+3241 


33 


16 


0. 


.485; 


RXCJ1122.2+6712 


22 


8 


0. 


.364; 


RXCJ1204.4+0154 


12 


9 


0. 


750' 


RXCJ1223.1+1037 


4 


2 


0. 


500; 


RXCJ1324.1+1358 


6 


2 


0. 


333; 


RXCJ1440.6+0328 


15 


7 


0. 


.467; 


RXCJ1604.9+2355 


9 


7 


0. 


778' 


All Groups 


146 


78 


0. 


534; 



+0.094 
0.084 
■+0.278 
0.231 
0.10.1 
0.084 
+0.181 
0.144 
,+0.226 
0.177 
+0.123 
'-0.091 
+0.051 
0.047 

+0.237 
0.173 
,+0.269 
-0.184 
,+0.231 
0.239 
+0.154 
-0.120 
+0.179 
-0.126 
,+0.250 
■0.245 
0.500 
0.323 
+0.440 
■0.215 
0.251 
-0.172 
+0.222 
-0.287 
-0.068 
-0.060 



109 


109 


80 


49 


41 


0.837; 


19 


75 


35 


6 


5 


0.833; 


67 


67 


26 


25 


15 


0.600; 


30 


30 


19 


15 


12 


0.800; 


23 


23 


17 


16 


14 


0.875; 


36 


36 


11 


18 


6 


0.333; 


284 


340 


188 


129 


93 


0.721 



20 
22 
13 
33 
22 
12 
1 
6 
15 
9 

156 



28 
22 
13 
33 
22 
12 
4 
6 
15 
9 
164 



16 
11 
10 
16 

8 
9 
2 
2 
7 
7 
88 



11 

7 
9 
9 
8 
7 
2 
3 
9 
3 
68 



10 

5 
8 
6 
5 
6 
1 
2 
5 
3 
51 



+0.152 
-0.130 
+0.167 
-0.360 
,+0.198 
-0.153 
+0.200 
-0.228 



-0.132 
+0.083 
-0.075 

0.909±°;°! 3 
0.714+pi 
0.889_ ; 308 

Rfi7+0.333 
u.oo<_ 264 

U.DZO_ Q 270 

0.8571^ 
0.500±°;£° 

U ' DD '-0.431 

556+ - 376 
UOOD _0.240 

1 000+ 000 
r.uuu_ 544 

u. 'OU_ 0il05 



53 


53 


44 


10 


10 


33 


28 


28 


17 


15 


15 


12 


16 


16 


14 


18 


18 


6 


170 


170 


126 


15 


15 


11 


12 


12 


9 


9 


9 


8 


9 


9 


6 


8 


8 


5 


7 


7 


6 


2 


2 


1 


3 


3 


2 


9 


9 


5 


3 


3 


3 


77 


77 


59 



Note. — Morphological and demographic information for our two samples. Columns arc: (1) Cluster or group name; (2) Number of 
objects successfully fit by GALFIT with Mr < — 20 mag; (3) The number objects from column 2 with n > 2.5; (4) Fraction of objects from 
column 2 with n > 2.5; (5) Number of confirmed objects with Mr < —20 mag; (6) Number of objects after a completeness correction, if 
any; (7) Number of objects in the sample with n > 2.5, corrected for completeness, if applicable; (8-13) Same as columns 2 -7, but with t he 
brighter magnitude cut (Mr — A-/£, + 1). The errorbars on the fractions are all single-sided, onc-sigma confidence intervals ( Gchrcls 1986). 
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Fig. 5. — Bulge-to-total flux ratio versus Sersic index for all of 
the groups and clusters in our sample. Points are as in Figurel2lfor 
X-ray AGN (filled circles), BPT AGN (stars), and inactive galaxies 
(open circles). This plot includes all the galaxies successfully fit 
with GALFIT, though the BPT AGN are only selected from the 
subsample with adequate Sloan spectra. The dashed line is drawn 
at n = 2.5. 




fracDev=0.66 



0.2 0.4 0.6 0.8 1 

SDSS fracDev 

Fig. 6. — Sersic index versus SDSS fracDeV parameter for all of 
the galaxies with SDSS coverage. Points are as in Figure [2] for X- 
ray AGN (filled circles), BPT AGN (stars), and inactive galaxies 
(open circles). The solid line is a linear fit and the dashed lines are 
drawn at n = 2.5, where we separate early and late-type galaxies, 
and the interception of the fit line with n = 2.5, at f racDeV = 0.66. 

for clarity. In both figures the top panel shows the results 
for the absolute magnitude threshold of Mr < — 20 mag 
and the bottom panel for Mr < M R + 1. These figures 
indicate that the AGN fraction is smaller in environments 
characterized by a higher velocity dispersion. 

Due to the small number of AGN in individual groups 
and clusters, the statistical significance of these trends is 
difficult to discern. We thus bin the data for all groups 
and all clusters separately, where we divide the two sam- 
ples at a velocity dispersion of a = 500 km s -1 , and 
compare these two environments. The binned results are 



presented in Table [7] and the right-hand panels of Fig- 
ures [7] and [H These results make the trend with velocity 
dispersion more clear. For the higher-luminosity host 
galaxies (Mr < M R + 1), the errorbars on the binned 
AGN fractions for the groups and clusters do not overlap 
for both all AGN and just those with early-type hosts. 
Specifically, for all galaxies more luminous than M R + 1 
we find that the X-ray AGN fraction is f A = 0.091±°;°g 

for groups and 0.047lg;Qig for clusters, or a factor of 
two higher in groups. The trend is somewhat more pro- 
nounced for the early- type galaxies, where the AGN frac- 
tion is f A ,n>2.5 = 0.119lg;gJ in groups and 0.048±g;g?| 
for clusters. This demonstrates that the AGN fraction 
is a factor of two higher in groups relative to clusters 
and that the AGN fraction is similarly higher when just 
early- type host galaxies are considered. 

The errorbars on these binned results, shown in the 
righthand panels of Figures [7] and [8j show the same 
single-sided, lcr confidence intervals as the results for in- 
dividual groups and clusters. In principle, if these binned 
error bars exactly overlap, then each population is dis- 
tinct from the other with 84% confidence, and if not, 
then the confidence level can be obtained by expanding 
or contracting the confidence limits until they exactly 
overlap. However, there is additional uncertainty due to 
the choice of a = 500 km s _ to separate groups and 
clusters. While this choice was physically motivated and 
is not an unreasonable point to divide the sample, there 
are additional, physically meaningful values of the veloc- 
ity dispersion to separate groups and clusters. For exam- 
ple, we could have chosen to divide groups and clusters 
at er = 400 km s _1 instead of 500 km s _1 . Therefore, 
proper statistical analysis of the difference between the 
two samples needs to include an additional penalty that 
reflects the other options for binning the data. We ac- 
count for this by raising the previous probability to an 
exponent that represents all of the other ways we could 
have binned the data. Thus, the probability that two 
populations are distinct from one another is 

(1 - (1 - CL) 2 ))^ ""** 8 , ( 3 ) 

where CL is the confidence limit at which the error bars 
just overlap, and A group i ngs is the number of possible 
groupings. For example, a confidence limit of 84% (one- 
sigma) and three possible groupings yields a 93% proba- 
bility that the two p opulations are distinct. Note that in 
iSivakoff et al.l (|2008l ) a similar analysis was performed 
to show that groups and clusters were different from 
one another, but there all possible combinations of the 
data were considered. In our case there are in principle 
(2 s ) = 120 possible ways to create two samples out of 
our data, but we do not adopt this approach as most of 
these options are nonconsecutive and not physically mo- 
tivated. From this analysis we conclude that the X-ray 
AGN and early-type X-ray AGN fractions are higher in 
groups relative to clusters at the 85% and 92% level, re- 
spectively, for the Mr < M R + 1 sample. The statistical 
significance for the Mr < —20 mag sample is lower (77% 
and 76%, respectively), although the trend is consistent. 

Another intriguing question we can begin to answer 
is how the AGN fractions in groups and clusters com- 
pare to the field value. While there is little data on 
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TABLE 7 
X-ray AGN Fractions 



Cluster / Group 


Nxjit,n>2.5 N x jit 


fx 


fx,n>2.5 


"X,fit,n>2.& X.fit 


fx 


fx,n>2.5 


(1) 


(2) (3) 


(4) 


(5) 


(6) (7) 


(8) 


(9) 



A85 

A644 

A3128 

RXCJ0110.0+1358 
RXCJ0746.6+3100 
RXCJ1022.0+3830 
All Clusters 



2 
2 
4 

1 
3 
12 



n nis+ - 024 
U.U18_ 012 

027+ - 035 

u - uz '-0.017 

n 060+ - 047 

U.UDU_q 02 g 

n ooo+ 061 

u.uuu_ 000 
0.043l°;J°« 
0.0831°;^ 
0.035±g;glg 



n 025+ - 033 
029+ 066 

U.U^»_ Q Q24 

u - llo -0 .063 

n nnn+ 097 

u.uuu_ 000 
0.059l°; 3 | 
0.182l°? 4 ° 



2 

1 

1 
2 
6 



0.0381°;°" 
0-0251°;^ 

U.UCSD_q 030 

o.oool ,;^ 

n 062+° 144 

U.UDZ_g 052 
U-J-D' -0.091 
u ' u4 ' -0.016 



U.U40_ Q Q2g 

000+ 056 
u.uuu_ 000 

0.059±g-.S|g 

o.oool ,;" 3 

u.uii_ 059 

n ooQ+0.440 
u.ooa_ 215 

0.048_ aol9 



A3125 
A89B 

RXCJ0844.9+4258 
RXCJ1002.6+3241 
RXCJ1122. 2+6712 
RXCJ1204.4+0154 
RXCJ1223. 1+1037 
RXCJ1324. 1+1358 
RXCJ1440.6+0328 
RXCJ1604.9+2355 
All Groups 



2 
2 
1 

2 





7 



107+ - 104 
u - iu ' -0.058 

091+° 120 
u.uai_ 05g 

077+ - 177 
u ' u( ' -0.064 

000+ 056 
u.uuu_ 000 

136+ 0133 
lou -0 .074 

n ooo+ 0153 
u.uuu_ 000 

o ooo+ - 460 

u.uuu_ 000 

000+ ' 307 
u.uuu_ 000 

0.000+ 123 



-0.000 
n+0.205 
-0.000 

n n55+ ' 025 

u.uoo_ 0i01g 



o.oocr 



n 195+ 0165 

u.izo_ 081 

0.182t° 24 ? 

100+ - 230 
u.iuu_ 083 

n ooo+ - 115 
u.uuu_ 000 

U - ZOU -0.161 

o 000+ 0205 
u.uuu_ 000 
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Note. — X-ray AGN fractions for the cluster and group samples. Columns arc: (1) Cluster or group name; (2) Number of X-ray 
identified AGN with good fits, an early-type morphology (n > 2.5), and Mr < —20 mag; (3) Number of all X-ray identified AGN with 
good fits in the sample; (4) Fraction of the fit galaxies that are X-ray identified AGN; (5) Fraction of the fit galaxies that are X-ray 
AGN and have early-type morphologies (n > 2.5); (6-9) Same as (2-5) but with the brighter absolute magnitude cut of Mr < M R + 1. 



X- ray selected AGN f ractions in the field, the study 
of iLehmer et al.l ()2006D examined the X-ray fraction in 
early- type galaxies as a function of redshift. When cal- 
culated with the same selection criteria as we employ, 
the field early-type AGN fraction is 6.6% for Mr < 
—20 mag (B. Lehmer, private communication, see also 
IMartini et all 120071) . Intri guingly, the group and clus- 
ter early-type AGN fractions are both consistent within 
the errorbars with the field value, although the cluster 
fraction is lower and the group fraction is higher. 

4.2. Emission-line and X-ray AGN 

The eight groups and three poor clusters selected from 
the NORAS sample have SDSS spectra as well as imaging 
data. Consequently, we were also able to classify them as 
BPT AGN or not based on their visible- wavelength emis- 
sion line ratios (see 5 l2~!T2"|) . We calculate the BPT AGN 
fraction as the ratio of BPT AGN in galaxies more lumi- 
nous than Mr = — 20 mag divided by the total number 
of galaxies above this luminosity with SPS S spectra. As 
noted above, there are 14 AGN above the IKewlev et al.l 
(2001) classification line of 199 total galaxies with spec- 
tra. Six of these BPT AGN are in 88 cluster galaxies and 
eight are in 111 group galaxies. The corresponding BPT 
AGN fractions are 0.068i^2 7 for groups and 0.072^°;°^ 
for clusters. These fractions are consistent. Only one of 
these 14 BPT AGN is also classified as an X-ray AGN 
and it is a member of a group (RXCJ0844. 9+4258). For 
comparison the X-ray AGN fraction is 0.039to'oi4, al- 
though is drawn from a larger host galaxy population 
than the two BPT fractions quoted above. Our BPT 
AGN fraction is comparable but somewhat lower than 
the X-ray AGN fraction of Ja ~ 0.07 (based on one ob- 
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Fig. 7. — X-ray AGN fractions versus velocity dispersion. The 
top panel of this figure shows this relationship for an Mr < —20 
cut o n host galaxy absolu te magnitude, as in Martini et al. (120061) 
and JSivakoff et al. 20 Q3). The bottom panel uses the c utoff of 
M R = Mjj + 1, as in lMartini. Sivakoff. fc Mulchaevl H2009I ). Mea- 
surements for all groups with a < 500 km s~ 1 (filled squares) and 
clu sters ( filled tria ngles) with at least one AGN are shown. Follow- 
ing Sivakoff et al. (2008), the right panels show the average AGN 
fraction for a < 500 km s _1 (squares), a > 500 km s _1 (triangles). 

ject) found by IShen et all (|2007D for M R < -20 mag. 
Their BPT AGN fraction is also ~ 7%, although this was 
calculated for Mr < — 18 mag. At this lower threshold 
they measure an X-ray AGN fraction of 0.7% (one out 
of ~ 140 galaxies). We similarly find BPT AGN that are 
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Fig. 8. — As in Figure [7] for the early type X-ray AGN frac- 
tions versus velocity dispersion. The top panel of this figure shows 
a measurement of t he field early-type X-ray AGN fraction from 
ILehmer et"aLl i200B ). 

not classified as X-ray AGN and X-ray AGN that are not 
classified as BPT AGN. 

We also compare the morphologies of the BPT AGN 
and the X-ray AGN to determine if the lower frac- 
tion of hi gh-luminosity AGN in denser regions found 
in SDSS (|Kauffmann et aLll2004f ) is correlated with the 
morphology-density relation. To make this comparison 
we plot the cumulative fraction of X-ray and BPT AGN 
as a function of Sersic index in Figure^ While this figure 
may suggest that the X-ray AGN are more likely to be 
in early- type host galaxies, both a KS test and a Mann- 
Whitney [/—Test indicate that the samples are consis- 
tent. The BPT AGN morphologies are also in very good 
agreement with the galaxies not classified as AGN by ei- 
ther method. The X-ray and BPT AGN shown in this 
figure are only the subset with spectroscopy from SDSS, 
and therefore not all of the X-ray AGN are shown. 

5. DISCUSSION 

The results of the previous sections have shown that 
there is tentative evidence that the AGN fraction is lower 
in clusters than in groups, and that this difference also 
holds when only early-type galaxies are considered sep- 
arately. This second statement is important because 
it helps to break the degeneracy between morphology 
and density and between morphology and AGN frac- 
tion. Thus the larger fraction of AGN in groups indi- 
cates that the AGN fraction is higher in both early and 
late-type galaxies and is not simply due to a larger frac- 
tion of late-type galaxies in lower velocity dispersion en- 
vironments. These results have interesting implications 
for AGN triggering and fueling, as well as the evolution 
of galaxies in groups and clusters. If true, one of the 
main implications is that early-type galaxies in groups 
are more likely to be AGN than their counterparts in 
clusters. This could be explained by higher cold gas frac- 
tions or greater likelihood of triggering, such as due to 
an interaction; most likely both play a role. In clusters, 
the higher rates of galaxy interactions and gas stripping 
could remove potential fuel from the galaxies more effi- 




Sersic Index 

Fig. 9. — Cumulative fraction of galaxies with Sersic index n 
less than the amount shown for the sample selected from the NO- 
RAS catalog. The BPT AGN sample (solid histogram) includes 
a larger fraction with small Sersic index, corresponding to later- 
type morphologies, while the X-ray AGN (dashed histogram) hosts 
are predominantly early-type galaxies; however, the difference is 
not statistically significant. All other galaxies with spectra and 
morphological fits are also shown (dotted histogram). There are 
14 BPT AGN and six X-ray AGN in these histograms. The last 
two X-ray AGN are at n = 5 and one AGN is common to both 
distributions. 



ciently than those same interactions driving angular mo- 
mentu m loss and AGN fue ling. It has also been pro- 
posed (|Hopkins et al. 2008a, bj) that smaller, less mas- 
sive groups are the ideal environment for AGN activ- 
ity caused by mergers. Though these studies focus on 
higher-luminosity AGN and rely on a different triggering 
mechanism than in this work, it would be interesting to 
investigate whether our weakly observed trend continues 
to even lower-density group environments discussed by 
iHopkins et all (|2008a| ). 

These results also help to explain the range of con- 
clusions that have been drawn about the AGN frac- 
tion as a fu n ction of environment. For example, 
iDressler et ail ()1985f ) found the BPT AGN fraction in 
clusters is a factor of five times lower than in the 
field. Many studies have confirmed this result for 
higher-luminosity AG N, yet found less of a difference 
at lo wer luminosities (IShimada et all l2000t [Miller et"aTI 
l200l iKauffmann et all 12004 iGrogin et alj I2005T) . In 
contrast, X-ray observations of clust ers find many more 
AGN than via the BPT technique (|Martini et all 120061 . 
l2007f K while observations of optically-selected, poor 
groups of galaxies identi fy a larger fractio n of BPT AGN 
relative to X-ray AGN (jShen et al.ll2007f l. Our spectro- 
scopic analysis of the rich groups and poor clusters se- 
lected from the NORAS catalog demonstrates that much 
of these pronounced differences in the literature are due 
to differences in selection because the BPT AGN and 
X-ray AGN are nearly disjoint populations. There are 
substantial populations of both types of AGN in groups, 
while this is less often the case in higher galaxy den- 
sity environments . The very low density groups in the 
iShen et all (|2007f ) sample are unlikely to be virialized 
systems, even though they are typical of groups found 
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via redshift surveys such as SDSS. The X-ray emission 
from the more massive groups and poor clusters sug- 
gests that these are virialized systems, and the change 
between mostly BPT AGN and X-ray AGN may reflect 
a change in the dominant accretion mode between unviri- 
alized and virialized sy stems. This is furthe r supported 
by the cluster sample of lMartini et alj (|2006f ). who found 
few BPT AGN and none that were not also detected as 
X-ray AGN, although many of these spectra were also 
of low signal-to-noise. For comparison, we find approxi- 
mately equal numbers of BPT AGN in groups and clus- 
ters, including examples of BPT AGN in clusters that 
are not X-ray AGN. Nevertheless, from the several pre- 
vious studies (jDressler et all Il999t iMartini et all [2006) 
the trend is that while the AGN population is lower in 
clusters than in groups, the decrease is larger for BPT 
AGN than for X-ray AGN. The more constant X-ray 
AGN fraction with local density is also sim ilar to the 
trend seen in radio AGN bv lBest et all (f2005h . who found 
that the fraction of radio AGN is relatively insensitive to 
environment, compared to BPT AGN. 

A potential physical explanation of these trends is 
that the BPT and X-ray AGN trace different accretion 
modes. The BPT AGN exhibit line ratios characteris- 
tic of AGN whose spectral energy distributions are well- 
matched by thin disk models where the accretion rate is 
greater than ~ 1% of the Eddington rate. In contrast, 
the weak or absent emission lines, at least in these mod- 
erate signal-to-noise spectra, combined with their sub- 
stantial X-ray luminosities, suggest radiatively-inefficient 
accr etion with lower accretion rates relative to Edding- 
ton dNaravan. Mahadevan. fc Quataertj|1998t [Ho1ll999l: 
iVasudevan fc Fabianl 120071 ). This simple picture is sup- 
ported by the weak trend that the X-ray AGN are more 
often found in early-type hosts than the BPT AGN. As 
both populations have comparable total luminosities, the 
X-ray AGN hosts have larger spheroids and are expected 
to have larger black hole masses as a result. As radio and 
X-ray emission are reasonably well correlated in the se ra- 
diatively inefficient models fe.g. lMerloni et afll2003f ). this 
hypothesis is also consistent with the distribution of the 
observed radio AGN fraction. 

Another interesting implication of this work is that 
the average AGN fraction for the groups and clusters to- 
gether is similar to the field early-type X-ray AGN frac- 
tion of 6.6±| |% (B. Lehmer 2006, private communica- 
tion), although the groups alone have a higher fraction. 
One interpretation of this result is that while galaxies in 
the field typically have substantial supplies of cold gas, 
perhaps even more than found in group galaxies, this 
is offset by a lower rate of triggering due to the lower 
density, at least if interactions and mergers play an im- 
portant role. This implies that the cluster environment 
is too dense, the field environment is too sparse, and 
the group environment is "just right" to fuel and trigger 
X-ray AGN. Further observations are required to draw 
firm conclusions about the relative AGN fraction in the 
field, groups, and clusters, yet if the variation we find 
is confirmed it would strongly point to the importance 
of galaxy interactions (although not necessarily mergers) 
for fueling even lower-luminosity AGN. This would be 
surprising, as it conflic ts with current studies based on 
companion fractions (|Fuentes- Williams & Stockd 119881: 



lSchmittll200lHLi et al.l l2008l. although those studies are 
all based on BPT AGN, rather than X-ray AGN. 

6. CONCLUSIONS 

The distribution of AGN as a function of environment 
is a potentially valuable probe of the fueling and trigger- 
ing of AGN, as well as the connection between galaxy and 
black hole evolution. We have conducted a new survey 
of eleven rich groups and poor clusters selected from the 
NORAS catalog at 0.02 < z < 0.06 to measure the AGN 
fraction as a function of environment. Our group sam- 
ple, defined to be environments with a < 500 km s~\ 
co ntains eight new grou ps plus two previously published 
in Sivakoff et al. (2008) and thus represents a factor of 
five increase in sample size. The cluster sample contains 
three new clusters and three previously published, or a 
factor of two increase in sample size. 

We identify X-ray AGN in these groups and clusters 
with a combination of spectral fits and flux ratio argu- 
ments to demonstrate that the X-ray emission from each 
X-ray AGN is inconsistent with other sources of X-ray 
emission. The main result of this analysis is that the 
X-ray AGN fraction is approximately a factor of two 
higher in groups than in clusters. This trend is ap- 
parent for both AGN in host galaxies more luminous 
than Mr — — 20 mag and for more luminous hosts with 
Mr < M R + 1, although the difference has greater sta- 
tistical significance for the higher luminosity threshold. 
The X-ray AGN fractions for the higher threshold are 
Ja{L x > 10 41 ; M R < M* + 1) = 0.047±g;gf| for clusters 

and 0.091±g;g|| for groups. There is a 85% probability 
that the group AGN fraction is larger than the cluster 
AGN fraction. This result may be more significant for 
the more luminous galaxy subsample due to the larger 
fraction of the most luminous galaxies that are X-ray 
AGN. 

Because the incidence of AGN in galaxies depends on 
host galaxy morphology, and the distribution of galaxy 
morphologies depends on environment, we have con- 
ducted the first quantitative morphological analysis of 
the AGN fraction in dense environments with these six 
clusters and ten groups. The morphological data for 
every confirmed group and cluster ga laxy was obtained 
with the GALFIT software package bv lPeng et all (|2002ft 
and used to separate early- type and late- type galaxies. 
We then calculated the X-ray AGN fraction for just the 
early-type galaxy populations in the groups and clus- 
ters separately and found that the early-type AGN frac- 
tion is also a factor of two higher in groups relative to 
clusters. For the higher galaxy luminosity threshold we 
find f A , n >2. 5 {Lx > 10 41 ;M fl < M* + 1) = 0.048±g;gf| 
for clusters and 0.119^044 for groups (92% confidence). 
The similar trends for early-type galaxies and all galax- 
ies indicate that the AGN fraction is not different simply 
because the morphological mix of galaxies changes as a 
function of environment, but rather that all galaxy types 
have a higher probability of hosting an AGN in the group 
environment. In addition, the group value is also higher 
than the best estimate of the early-type AGN fraction 
in the field. This may be because group galaxies, even 
early- type group galaxies, are more likely to have sub- 
stantial cold gas reservoirs for AGN fueling than cluster 
galaxies, while galaxy interactions are more likely to oc- 
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cur in groups than the field, or some combination of these 
effects. 

Finally, we have also estimated the AGN fraction based 
on emission-line diagnostics for the subset of the galax- 
ies with SDSS spectroscopy. There are 14 BPT AGN 
in this subset of groups and clusters, as well as six of 
our nine X-ray AGN. Strikingly, these populations are 
nearly completely disjoint: only one AGN meets our cri- 
terion as both an X-ray AGN and as a BPT AGN. This 
is a clear demonstration of how different selection tech- 
niques may identify different populations of AGN. Our 
morphological analysis shows that the host galaxies of 
these two AGN types are marginally different in that the 
X-ray AGN are more likely to be hosted by early-type 
galaxies. While the host galaxies for both AGN pop- 



ulations are more luminous than Mr < —20 mag, the 
earlier-type hosts for the X-ray AGN imply relatively 
larger supermassive black holes compared to the BPT 
AGN. These observations are thus consistent with lower- 
efficiency, but relatively more X-ray bright, accretion in 
the X-ray AGN. 
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